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ABSTRACT 

We present the first dynamical mass measurements for Lyman-a galaxies at 
high redshift, based on velocity dispersion measurements from rest-frame optical 
emission lines and size measurements from HST imaging, for a sample of nine 
galaxies drawn from four surveys. These measurements enable us to study the 
nature of Lyman-a galaxies in the context of galaxy scaling relations. The result- 
ing dynamical masses range from 10 9 to 1O 1O M . We also fit stellar population 
models to our sample, and use them to plot the Lyman-a sample on a stellar 
mass vs. line width relation. Overall, the Lyman-a galaxies follow well the scal- 
ing relation established by observing star forming galaxies at lower redshift (and 
without regard for Lyman-a emission), though in ~ 1/3 of the Lyman-a galax- 
ies, lower-mass fits are also acceptable. In all cases, the dynamical masses agree 
with established stellar mass-linewidth relation. Using the dynamical masses as 
an upper limit on gas mass, we show that Lyman-a galaxies resemble starbursts 
(rather than "normal" galaxies) in the relation between gas mass surface density 
and star formation activity, in spite of relatively modest star formation rates. 
Finally, we examine the mass densities of these galaxies, and show that their 
future evolution likely requires dissipational ("wet") merging. In short, we find 
that Lyman-a galaxies are low mass cousins of larger starbursts. 

Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: 
evolution 
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Introduction 



Lyman a line emission is an increasingly important tool for identifying actively star 
forming galaxies in the distant univ erse. First proposed as a potential signpost for primitive 
galaxies (IPartridge fc Peebles! 1 19671 ). Lyman-a emission has now been used to identify thou- 
sands of galaxies at redshifts 2 < z < 7, along with smaller samples at z < 2 and several 
candidates at z > 7. 

While no galaxy yet identified by any means is demonstrably primordial, Lyman-o; 
selected samples do have several properties suggestive of youth. Thei r starlight is domi- 
nated by young popul ations with chara cteristically low stellar masses (IPirzkal et al.l 120071 : 
Finkelstein et all l2007h and small sizes faond et all bood : iMalhotra etaPboij l. Yet, the 
correlation properties of these object s sugge s t that they are assoc iated with moderately large 



halos (mass ~ 10 M ; iKovac et al.l (120071 ); iGuaita et al.l (120101 ) ) . Combining these results 



suggests that a Lyman-a galaxy contains only a small fraction of the baryons that should 
be associated with its host dark matter halo. It would be interesting to know whether the 
"missing" baryons are present, either as old stellar populations, or in the interstellar medium 
of the Lyman-a galaxy. 

Dynamical mass estimates for Lyman-a galaxies could potentially address this question, 
providing a standard for comparison with both the stellar masses and the dark halo masses. 
However, such dynamical mass estimates require an accurate measurement of the galaxy's 
velocity dispersion. The easiest approach to kinematics would be to use the Lyman-a line 
width, which is measured for most spectroscopically confirmed Lyman-a galaxies. Unfortu- 
nately, this does not lead to useful velocity dispersion information: The Lyman-a line profile 
can be dramatically affected by the interplay of resonant scattering and gas kinematics in 
the emitting galaxy. 

We therefore turn in this paper to studying Lyman-a galaxy kinematics using the strong 
rest-frame optical emission lines of [OIII] A5007A and Ha A656lA. We buil d on the first de- 
tections of such lines in Lyman-a selected galaxies at redshifts 2.2 < z < 3.1 ( jMcLinden et al. 
201 ll ; iFinkelstein et al.ll2011t lHashimoto et al.ll2012f ). In section |2l we describe the sample 
and the key data. In section [3] we estimate dynamical masses based on the observed line 
widths. In section HI we analyze the stellar populations and dust reddening in these galax- 
ies, using deep archival photometry. In section we combine the rest-optical line width 
measurements with the stellar masses from population synthesis modelling to compare these 
Lyman-a galaxies with expectations from the stellar mass Tully-Fisher relation. In section |6l 
we examine the relation between gas mass surface density and star formation surface den- 
sity to show Lyman-a galaxies lie on the same sequence as starburst galaxies. Finally, in 
section we explore the mass densities of the sample and the implications for the future 
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evolution of Lyman-a galaxies. 

Throughout the paper, we adopt a A-CDM "concordance cosmology" with fl M = 0.27, 
Q A = 0.73, and H = 71 kms" 1 Mpc" 1 . 



2. Description of Sample 

Our sample consists of Lyman-a emitting galaxies that were selected from four surveys. 
All are selected directly by the presence of a strong Lyman-a emission line in the survey 
data. For the analysis in this paper, we make use of (a) velocity dispersions a, derived 
from spectra of rest frame optical emission lines; (b) sizes, as measured by the half-light 
radius r e in broad-band optical Hubble Space Telescope images; (c) stellar masses and dust 
extinctions, derived from spectral energy distribution (SED) fits; and (d) star formation 
rates (SFR), determined from Ha line flux measurements where available (and from spectral 
energy distribution fitting otherwise), and corrected for dust using the results of the SED 
fits. In this section we summarize the sources for various galaxies in the sample, along with 
the measurements of a, r e , and SFR. The stellar mass derivations are discussed later, in 
section @J The properties of the sample are also summarized in tables 1-3. 



2.1. Bok telescope z — 3.1 survey objects: 



First, we select two Lyman-a galaxies from a 5020A narrowband survey using the Stew- 
ard Observatory's 90 inch Bok Telescope on Kitt Peak, Arizona. The objects we study here 
were spectroscopically confirmed as Lyman-a emitters at z ~ 3.12 using the 6 .5m MMT 
on Mt. Hopkins, Arizona. This survey was introduced in McLinden et al. (boil), an d fur- 
ther details will be presented in a forthcoming paper (IMcLinden fc all 12013. in prepl ). We 
have followed up galaxies from this sampl e with three near-inf rared spectrographs: LU- 
CIFER, on the Large Binocular Telescope (IMcLinden et al.ll201ll ); NIFS, o n Gemini North 
(|Richardson et al.l 12013. in prepl ); and NIRSPEC, on the Keck II telescope (IMcLinden fc al 
2013. in prepl). The objects stud ied here are LAE40844 and LAE27878 (following the num- 
bering in IMcLinden et al.l (120111 )). 



The prim ary source of line widths for these objects is our Gemini NIFS data (IRichardson et al 
2013. in prepl ). since our Keck observations of 40844 were primarily aimed at the (still un- 
detected) [Oil] AA3726,3729 lines, while the kinematic line width is not well resolved in the 
LUCIFER observations. These Gemini observations yielded line widths of AVfwhm = 13. 8A 
for object 40844 and 8.6A for object 27878. The instrumental resolution was about 5A. Sub- 
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tracting this instrumental resolution in quadrature from the measured line widths yields 
13. OA and 7.0A respectively. The corresponding line-of-sight velocity dispersions become 
crj os = 80km s" 1 and crj og = 43km s -1 respectively. 

The half-light rad ii for these objects are 1.1 kpc for LAE40844 and 1.3 kpc for LAE27878 
(IMalhotra et al.ll2012l ). We estimated their star formation rates from SED fitting (section HJ), 
since their Ha line falls at 2.70/xm and cannot be easily observed. We obtained 12OM yr _1 
for LAE40844, and 34M yr _1 for LAE27878. We estimate the uncertainties in these star 
formation rates from the range of SFR in models yielding acceptable fits to the data (see 
section HJ) . 



2.2. HETDEX Pilot Survey objects: 



Second, we use two Lyman-a galaxies selected using a blind spectroscopic search with 
an integral field spectrograph as part of the HETDEX Pilot Survey. These are objects 
HPS194 (z = 2.287) and HPS256 (z = 2.491). For these, we base our kin e matic line 
widths on Kec k+NIRSPEC Ha line measurements from iFinkelstein et al.l (120111 ) ; ISong fc al 
( 2013. in prepl ). These used low-resolution mode with the NSPEC-7 blocking filter. The 
measured line widths were 18.1 ±0.6A (FWHM) for HPS 194, and 19.4±0.9A (FWHM) for 
HPS 256. The instrumental resolving power is R ~ 1500 with the 0.76" slit, corresponding to 
14-15A FWHM, and direct measurements of sky lines in the two spectra yield resolutions of 
14.4A and 14. lA (FWHM) respectively. Subtracting these in quadrature yields line widths of 
11. OA and 13. 3A, respectively. The corresponding line-of-sight velocity dispersions become 
^los = 65km s -1 and crj og = 74 km s -1 respectively. T he half-light radii, as reported in the 



COSMOS ACS i-band catalog (ILeauthaud et al.ll2007f ). are 1.5kpc for HPS194, and l.lkpc 
for HPS256. (HPS 194 in fact corresponds to a pair of continuum sources separated by 
~ 0.5" in the COSMOS survey's HST+ACS images. We assume that the strong Lyman-a 
and Ha emission come from the brighter and more compact source. The other source has 
a 2.4kpc half-light radius, which would raise our dynamical mass for this object by about 
60%.) 



The Ha line fluxes (ISong fc all 12013. in prepl ) yield star formation rates of 17M Q yr 
and 2OM yr _1 , respectively, before dust correction. This is based on the conversion SFR = 
4.6 x 10~ 42 (Lff n./ ergs'^M^ .yr -1 , which is appropriate for a Chabrier (2003) initial mass 
function [IMF] ( ITwite et al.ll2012l ). The inferred SFR would be 1.2 x greater for a Kroupa 
(2003) IMF, and 1.8x greater for the IMF assumed in Kennicutt (1998). We correct these 
star formation rates for dust extinction within the emitting galaxy. The uncertainty in the 
star formation rate is dominated by the uncertainty in this dust correction, which we estimate 
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by considering the full range of extinction among models with \ 2 — Xmin + 1' anc ^ further 
allowing the possibility that the extinction of Ha could exceed the extinction of continuum 
starlight by up to a factor of 2. 



2.3. Subaru NB survey objects: 



Third, we use four Lyman-q ga laxies from a narrowband survey described by lHashimoto et al 

(120 121 ) and iNakajima et al.l (120121 ) . These objects are drawn from two fields — COSMOS 



(objects COSMOS13636 and COSMOS30679) and the Chandra Deep Field South (objects 
CDFS3865 and CDFS6482). We obtained line widths for all so urces by measuring the plot- 
ted FWHM of emission line profiles in figures 1, 2, and 5 of lHashimoto et al.l ( 120121 ). In 
general, the lines in this sample appear marginally resolved. In our analysis we regard the 
line width measurements as upper limits where appropriate. 

The COSMOS field spectra were from Keck + NIRSPEC spectroscopy. Here we ob- 
tained directly measured FWHM of 16.6 and 16.5A for COSMOS13636 and COSMOS30679, 
respectively. With a resolving power of R = 1500, the instrumental resolution corresponds to 
13. 8A, which we subtract in quadrature to yield nominal velocity dispersions o = 57km s _1 
and a = 55km s _1 , respectively. In each case, the measurement remains consistent with a 
fairly broad range, < a < 96km s _1 . 

For these COSMOS field sources, we base the half light radii on the COSMOS ACS 
i-band catalog, obtaining 0.79kpc and 1.82kpc respectively. Finally, their star formation 
rates as inferred from their Ha fluxes are 8.7M yr _1 and lO.5M yr _1 respectively. 



The CDFS field source spectra (lHashimoto et al.l 120121 ) were obtained with Magellan 
+ MMIRS, and have a somewhat lower resolving power (R « 1120, corresponding to 18A). 
The measured FWHM are 23A for CDFS3865 and 20A for CDFS6482. Subtracting the 
instrumental resolution in quadrature yields velocity dispersions of 105km s -1 and 42km s _1 , 
respectively, where again there is considerable uncertainty for the narrower, semi-resolved 
line (consistent with < o < 71km s _1 ). The star formation rates for these sources, based 
on their Ha line fluxes, are 125M yr _1 and 31M Q yr _1 respectively. 

For the physi cal sizes of the C DFS sources, we downloaded archival HST imaging from 
the GEMS survey faix et al.ll2004l) and measured the half light radii. We used the SExtractor 
faertin fc Arnoutsl Il996h half light radius, and as a consistency check also measured the 
fluxes in a series of circular apertures and interpolated the resulting photometric growth 
curve. Both methods gave consistent answers, with half light radii of 0.96kpc for CDFS3865, 
and 1.78kpc for CDFS6482. 
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2.4. ESO z = 2.25 survey object: 



Finally, we use one Lyman-a galaxy, L AE-COSMOS - 47, fr om a narrowband-selected 
z = 2.25 COSMOS field sample obtained bv iNilsson eTaD f|201lh using the ESO 2.2m tele- 
scope, with fol lowup observations obtained by Fynbo and collaborators using the X-Shooter 
spectrograph (jVernet et al.l 1201 ll ) on the VLT. The X-Shooter spectrum provided a well 
constrained velocity dispersion measurement of 30km s _1 . The star formation rate, derived 
from the Ha line flux, i s 33M yr _1 . The half-light radius, from the COSMOS ACS catalog 
( jLeauthaud et al.l 120071 ). is 1.14kpc. 



3. Dynamical Mass Estimates 

The measured velocity dispersions a from the rest-optical emission lines are a good 
estimate of the total luminosity-weighted kinematics of the gas. The galaxies we are studying 
are spatially unresolved even in ~ 0.5" seeing, meaning that the ground-based spectra we 
use effectively sample the integrated light, with no important dependence on slit width. 

The precise conversion from velocity width to mass will depend on the kinematic struc- 
ture of these galaxies. For a pure rotation-supported model with a flat rotation curve, we 
expect a 2 = sin 2 iv 2 /2, where i is the inclination angle of the disk (with i = 90° correspond- 
ing to an edge-on system). 

The simplest dynamical mass estimate from these measurements, which is a lower bound 
to the true gravitating mass, is Md yn > v 2 r/G, where r is the maximum radius at which we 
observe light from the galaxy. A more practical choice of radius is the effective radius r e , 
defined as the radius that encloses half of the galaxy's light in projection. If we presume that 
the half-light radius is also the half-mass radius, our revised estimate of the dynamical mass 
becomes M dyn rs 2v c (r e ) 2 r e /G ~ 4a 2 r e /(G sin 2 i) > Aa 2 r e /G. The resulting mass estimates 
range from 10 9 to 10 10 M Q , and are summarized in table [8) 

Several circumstances could affect our estimated mass. A rotating disk with i < 90° 
would reduce the measured a (relative to the edge-on case). Also, in our spatially unresolved 
spectroscopy, the observed a 2 reflects the luminosity-weighted average kinematics. If a sig- 
nificant fraction of the galaxy's light is emitted from regions where the local circular speed 
v c (r) is below the maximum circular speed (v max ), we should expect the weighted average 
o 2 to underestimate v max and hence the mass. The precise magnitude of this effect depends 
on the galaxy's light profile and rotation curve. Also, like any dynamical mass based on 
luminous tracers, our estimate is insensitive to mass located outside the luminous matter 
distribution of the galaxies. If the galaxies are embedded in extended dark matter halos, the 
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total mass of the halo could be many times the mass estimates derived from the observed r e 
and a 2 . 

Turbulence in the galaxy's gas would contribute to the measured a, though in virial 
equilibrium, that turbulence would constitute a source of pressure support and the ordered 
rotation of the galaxy would be correspondingly reduced. Most interestingly, if the galaxy 
is not in equilibrium at all, our assumed relations between kinematics and mass could be 
substantially wrong. Our mass estimate implicitly assumes that the virial theorem is fulfilled, 
that is, that the kinetic energy K and potential energy U of the galaxy are related by 
K = —U/2. On the other hand, in a cold accretion scenario, new material falling into a 
galaxy for the first time should have K = —U, i.e., the motions are faster for the same 
gravitating mass under these conditions, and the mass inferred from gas motions would be 
correspondingly over-estimated by a factor up to ~ 2. 

Given these uncertainties, it is best to regard our direct dynamical mass estimates as 
approximate numbers, good to a factor of perhaps 2 when regarded as lower bounds to the 
true dynamical mass. Other dynamical mass estimates in the liter a ture consider a more 



general scali ng coefficient so that M dyn = f3a 2 r P /G (see iToft et al.l (120121 ) and references 



therein, esp. Jorgensen et al. ( 1996 ) and Cappellari et al. ( 20061 )). These works favor (3 ~ 5 



for early-type galaxies with Sersic index n w 4, and find that despite theoretical expectations 
for some increase of with decreasing Sersic n, the observational evidence favors /3 ~ 5 for 
a wide range of n. Thus, the simple arguments that led us to use /3 = 4 likely come fairly 
close to the correct dynamical masses. 

A complementary approach to interpreting the kinematic data on these galaxies is to use 
their linewidths to place them on some form of the Tully-Fisher relation, and so to compare 
them on an equal footing to other galaxy populations. Such an approch avoids the difficulties 
associated with identifying the right radius to use in estimators of the form M ~ v 2 R/G. 
For high redshift galaxy populations, a small scatter with weak redshift evolution has been 
demonstrated for the stellar mass Tully-Fisher relation, and we place our galaxies on such a 
relation in section [5] below. To do so, we first need their stellar masses. 



4. Population Synthesis Modeling 



All of the galaxies we study have extensive multiband photometry in the literature, gen- 
erally including mu ltiband optical da ta, some deep ground-based near-infrared photometry, 
and Spitzer IRAC (IFazio et al.l 12004 ) observations that are deep enough to be constraining 
in at least the 3.6/xm channel. 
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We have used these data to derive stellar mass estimates for the full sample. Stel- 



ies ( 


Finkelstein et al. 


2011; 


2013. in prep 


). While these 



20121 ; iNakaiima et al.l 12012 ; McLinden fc al 



lap, we opted to fit the entire sample using a single procedure to avoid potential difficulties 
comparing masses derived using different methodologies. 



4.1. General comments on SED fitting: 

(1) The strong Lyman-a emission in these objects requires the presence of a young 
stellar population, whose ultraviolet light ionizes interstellar hydrogen that then recombines 
to produce the observed Lyman-a radiation. (2) The amount of dus t in the fitting is es - 



sentially determined by the UV spectral slope (cf. iMeurer et al.l (119971 ); lHathi et al.l ( 120081 ); 



Finkelstein et al.l (12011I )). (3) In some objects, no stars older than 10 7 years are required to 
explain the observed light. (4) A considerable mass in old stars is permitted. (5) The ionizing 
photon production for the best-fit stellar populations can be converted to a Lyman-a lumi- 
nosity by assuming that 2/3 of the ionizing photons are ultimately converted to Lyman-a 
radiation (as expected under Case B recombination with a negligible escape fraction for ion- 
izing photons). Combining this line luminosity estimate with the directly measured Lyman-a 
line flux gives an estimate of the escape fraction for the Lyman-a photons. These escape 
fractions are sensitive to the details of the star formation history over the last ~ 6 Myr. 
We can say with reasonable confidence that the resulting escape fraction is of order half for 
the most plausible models, and that Lyman-a escape fractions below 20-30% are ruled out 
unless we change the stellar population in some way that dramatically increases the ionizing 
photon production. (6) The estimated stellar population ages and masses would increase 
dramatically if we did not account for the [OIII]AA4959,5007 line fluxes. In this case, the 
stellar population fitting code attempts to interpret the red H— K color and the bright flux 
in the rest-frame 5000A range as due to older stellar populations. This effect can exceed an 
order of magnitude in both age and mass when the filter containing the [OIII]AA4959,5007 
lines is the reddest filter considered, while it is usually smaller when an additional filter 
redward of the 4000A break is included (meaning, in our case, the IRAC photometry). 



4.2. Starburst99 modeling: 



We measured stellar mas ses for our sample using the Starburst99 population synthe- 
sis code (ILeitherer et al.lll999[ ). For each galaxy we started with publically available broad 
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band photometry. For t he seven objects in the COSMOS field, w e used COSMOS project 
broad band photometry (ICapak et al.ll2007l ; iMcCracken et al.ll2010l ). either directly from the 



COSMOS archive, or as quoted in the paper s defining the samples (IFinkelstein et al.l 12011 



Hashimoto et al. 


2012 


Nakaiima et al.l 


2012) 


survey photometry ( 


Gawiser et alJ2006^ 


(aga 



20121 ) . For the two CDFS objects, w e used MUSYC 



( 120121 )). In all cases, our final model fitting used at least 9 photometric bands spanning at 



least the wavelength range OAfim < A < 3.6/im. 

To interpret the photometry in terms of stellar population parameters, we first ran 
Starburst99 to generate a grid of model spectra for stellar populations at a range of ages, 
from 2 x 10 6 years (which is so young that no star has yet left the main sequence), up to 
2 x 10 9 years (which is the age of the universe at z ~ 3.1), and using a Kroupa (2003) 
IMF. We then assumed a star formation history (as described below) and generated a model 
spectrum by a linear combination of the spectra for particular age steps from the Starburst99 
output. 

We accou nt for the mean opacity of the intergalactic medium using the prescription of 



Madaul ( 119951 ) . We do not treat the variance in the IGM opacity in the present work. (Doing 



so would effectively add uncertainty to the expected fluxes in the u*, Bj, and g bands, and 
so would improve the model x 2 ? but would not likely change the best fit model parameters 
much.) 

We model du st i n our sample galaxies using the Small Magellanic Cloud (SMC) ex- 
tinction law from iPeil ( 119921 ). and treating the extinction as a thin screen. The SMC law 
is a reasonable choice since its metallicity corresp onds approximately to the few metallicity 
constraints so far available for Lyman-a emitters ( IFinkelstein et al.l 1201 if ). 



We next added the directly measured spectroscopic line fluxes to the model spectrum at 
the appropriate wavelengths, since these emission lines are not included in the model output. 
This step follows the application of IGM and dust opacity, since nature has already applied 
these effects to the observed line fluxes. 

At this point, we have a full model spectrum accounting for stellar populations, emis- 
sion lines, dust, and intergalactic hydrogen absorption. We multiply this spectrum by 
the bandpass of each filter in the photometric data set, integrate, and normalize appro- 
priately to obtain model fluxes in each observed filter for the model under consideration. 
These can be compared to the observed data to obtain a goodness-of-fit parameter \ 2 — 

^2j=l(fj,obs fj,mod) I fj,obs) ■ 

To optimize the model, we fitted the observed spectral energy distribution of the sample 
galaxies by a simple Monte Carlo approach that randomly varies the amount of dust and 
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the mass in stars in each of 13 logarithmically spaced age bins, and accepts a change to 
the model parameters using x 2 minimization. This approach allows a more general star 
formation history than a single burst or an exponentially decaying star formation rate. While 
the resulting sampling of stellar ages is somewhat coarse compared to a typical single burst 
model, the associated uncertainty is not a dominant factor in our stellar mass estimates. The 
model fits usually converge to a case where 1-3 of the 13 mass bins dominate the luminosity 
at all observed wavebands. We therefore "trim" the parameter list by fixing the stellar mass 
to zero in bins that are clearly of minor importance, and rerun the fit with only the important 
mass bins. This "trimming" step generally has a negligible effect on the final x 2 i confirming 
that 1-3 simple stellar populations can explain the observed spectrum as well as a more 
complex star formation history. 

To further explore the parameter space of acceptable fits, we modified our code to 
optimize for either minima or maxima of either stellar mass or dust extinction, subject to 
the constraint that the model x 2 remain close to the \ 2 of the best fitting model for each 
object. We explored models with A% 2 = +1 and Ax 2 = +4. This corresponds approximately 
to the la and 2a error regions in the parameter space. 



5. The Stellar Mass Tully-Fisher Relation 



The stellar mass Tully-Fisher (SMTF) relation is a correlation between the kinematic 
line widths of galaxies and their stellar masses. The SMTF relation is more robust to 
differences in stellar population mass-to-ligh t ratio than the origi nal Tully-Fisher relation 
(which correlat es luminosity with l ine width; iTully fc Fisher! (11977| )). The relation is further 
generalized by Kassin et al.l ( 120071 ). who demonstrated that replacing the circular speed V c 

1 /2 

with the kinematic estimator S0.5 = (0.5V^ + a 2 ) results in an SMTF that is both tighter 
and more applicable to the wide range of galaxy properties seen at high redshift. So. 5 has an 
additional advantage: For spatially unresolved galaxies (like those we study here), S .5 can 
be measured reasonably accurately regardless of whether the width is dominated by ordered 
rotation or by random motions. 

Our actual measurement is a single number, the line width, characterized by the ob- 
served line-of-sight velocity dispersion a D b s . We then use S .5 — &obs for our galaxies. Consider 
pure circular motion with a flat rotation curve of velocity V c , viewed edge-on: we will find 
a obs = ((sin(0) \4) 2 ) = 0.5V 2 , so that So. 5 = a t, s . If instead the motion is entirely random, 
with V rot = 0, we expect a Q b s = a, and again, a a b s = So. 5. Only for ordered rotation in a 
face-on configuration do we expect a b s to be a significant under-estimate of S0.5. While we 
cannot rule out this possiblity with our data, it is likely that these galaxies are dynamically 
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hot (a/V ro t 1), since the stellar populations dominating the observed light are at most 
~ 1 dynamical time old. 

Our results are shown in figure (U both for stellar masses from SED fitting and for 
dynamical masses from line width and spatial extent. The error bars for the stellar masses 
show the ranges of stellar mass permitted by models with Ax 2 = X 2 ~ Xmin < 4, and with 
Ax 2 < 1. The error bars for dynamical mass are determined by the uncer tainties in a, and 



are di agonal since the x-axis of the plot is a. For three objects from the lHashimoto et al. 



(120121 ) sample, a and M dyn are plotted as upper bounds. 



Comparing our Lyman-a sample to the stellar mass Tully-Fisher relation that lKassin et al. 



(120071 ) reported for 0.1 < z < 1.2 star-forming galaxies, we find that both the stellar masses 
and dynamical masses of LAEs generally follow the established relation. Dynamical mass 
estimates show less scatter than stellar masses, as might be expected given the vagaries of 
star formation histories. 

About half the galaxies are consistent with stellar masses falling a factor of two or more 
below their dynamical masses, when we account for the range of acceptable estimates for 
both. For these galaxies, it is possible that the dynamical mass within the central 1-2 kpc is 
dominated not by the young stars that power the observed Lyman-a emission, but by some 
other unseen component. Collisionless dark matter should not be so strongly concentrated 
in the central kpc of a dark matter halo. Old stars (formed from centrally concentrated 
gas 3> 10 8 years ago) are a possible alternative. The dynamical mass estimates we present 
are in fact a tighter limit on the total mass in old stars within the inner ~ kpc than are 
the photometric limits. Given the actively star-forming nature of these Lyman-a emitting 
galaxies, a reservoir of gas is the most intriguing possibility for the "excess" dynamical 
mass. Overall, though, such excess mass is merely permitted and not required by the data- 
so, overall, the Lyman-a galaxy sample shows consistency with the stellar mass Tully-Fisher 
relation for other samples. Whatever physical properties allow Lyman-a to escape these 
particular star forming galaxies, they do not strongly affect the mass-linewidth relation. 

The tight correlation observed between M^„ n and line wi dth is related to the s mall and 



nearly constant physical sizes of the Lyman-a galaxy sample (IMalhotra et alj|2012j ). A fixed 



physical size, combined with variable (and sometimes uncertain) line widths, can generate 



the observed slope of the Md yn - So. 5 relation. The M* - So. 5 relation from lKassin et al.l (120071 ) 
has a somewhat steeper slope. Presuming that stars form a fairly large and fairly constant 
fraction of the dynamical m a ss, th is slope could be interpreted as evidence for size-linewidth 



relation in the iKassin et al.l (120071 ) sample, with r e oc Sq 5 for 7 ~ 1. 
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Fig. 1. — Stellar and dynamical mass Tully-Fisher relation. Colored points mark the Lyman- 
a emitting galaxies for which we present dynamical mass estimates. For each, we plot the 
best fitting stellar mass (blue) and the dynamical mass (green). The vertical bar through 
each point-pair spans the range of stellar masses for acceptable models (thicker line: models 
with Ax 2 < 1; thinner line: models with 1 < A% 2 < 4). The error bars on dyamical 
masses are diagonal, since uncertainty in velocity dispersion affects both plotted quantities. 
Point shapes correspond to surveys: Circles O Bok t e lescop e z — 3.1 survey; pentagons -H- 
HETDEX pil o t surv ey; triangles <B- Hashimoto et al.l ( 120121 ) Subaru survey; and square 4-> 
Nilsson et al.l (120111 ) survey. Objects with marginally resolved lines are plotted using their 
best-estimate values of a, but can be regarded as upper limit s on both Sp.z and M dyn . Dashed 



lines are the best-fit stellar mass Tully-Fisher relation from lKassin et al.l (120071 ) . 
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6. Star formation scaling laws 

We now take the dynamical mass as an upper bound on the gas mass in these galaxies, 
and compare their properties to the scaling relations that describe star formation in other 
galaxy classes. 

We use the gas surface mass density limit H gtmax = M<^„/(27rrg). We can improve 
this bound by subtracting our stellar mass estimates. In practice, our maximum stellar 
mass always exceeds the dynamical mass, allowing the possiblity that there is no gas mass. 
In most cases, though, our minimum stellar masses are below the dynamical masses, and 
we can subtract them to yield refined estimates of the gas mass surface density (£ 3 < 
(M dyn - M^ mm )/{2nr 2 e )). 

The star formation rate and gas surface density can be related according to scaling 
laws of the form log ~ c^logS^ + (3sf (where T, S fr is in M M yr -1 kpc^ 2 , and S fl is 



i n M^ kpc 2 ). For nearby spirals, and for distant star-forming BzK galaxies, iDaddi et al. 



(120101 ) find asF = 1-42 and (3sf = —9.83 for our choice of units. For submillimeter galaxies 
and (U)LIRGS, Daddi et al find a parallel but offset "starburst" sequence, with (3sf ~ —8.93 
(corresponding to 8x more star formation for the same gas surface density). 

We determined T>sfr for our sample using the SFR estimates and half-light radius mea 



surements discussed in [2j We have selected a Chabrier IMF for consistency with IDaddi et al. 



( 120101 ) . These values of T'Sfr exceed the expectations for "normal" star forming galaxies by 
a median factor of 4, based on S 9imax alone (i.e. assuming that all the gravitating mass is 
gas). The disagreement is sign ificant at the > 3cr level, relative to the 0.33 dex scatter in 



the scaling relation reported by lDaddi et al.l (120 lOh . A factor of 4 would place the Lyman-a 



galaxies in between the normal and starburst sequences, though closer to the starburst se- 
quence. If we use S 9 = {Md yn — M+ tmin ) / (2iTr 2 ) in the scaling relations, we find that the 
median galaxy in our sample is forming stars at twice the rate expected even under the 
starburst scaling. The bottom line from this comparison is that Lyman-a galaxies likely 
belong to a family of starbursting objects that includes ULIRGS and submillimeter galaxies, 
despite order-of-magnitude differences in mass and star formation rates. 

While the "normal" and starburst galaxies follow distinct J^sfr- &g relations, they obey 
a single relation when E 5 is replaced with the quantity S 9 /r dj/n , where r dyn is the dynamical 
time. We use r dyn « 2^r e /a to place our Lyman-a galaxy sample on this relation also. 
While the Lyman-a galaxies appear less unusual whe n measured against this relation, they 



remain systematically above the trend line found by IDaddi et al.l ( 120101 ) . The difference is 
suggestive, rather than significant, being a 2a effect. The median offset is a factor of ~ 3 in 
T, S fr at fixed T, g /Td yn - This is comparable to the 0A4dex scatter in the relation as reported 
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Fig. 2. — Star formation law comparison. First panel: Relation between star formation rate 
surface density T, S fr and gas mass surface densit y E g , both fo r our samples (large colored 
points) and for comparison samples drawn from iDaddi et al.l ( 120101 ) . For each Lyman-oj 
galaxy, the green point marks the upper bound on gas surface density obtained by associating 
the entire dynamical mass with gas. Cyan and blue points take a portion of the dynamical 
mass to be associated with stars, using the "2cx" and "la" low-mass models from out stellar 
population fitting. Error bars in YIsfr account both for uncertainties in Ha flux and in 
extinction corrections. Point styles identify subsamples, as in figure [H The solid line marks 
the relation for normal star forming galaxies, and "x" points the galaxy samples obeying 
that relation. The dashed line marks the starburst relation, and open black circles mark 
the corresponding galaxies. The Lyman-o: galaxies are inconsistent with the normal star 
formation relation, and consisten t with the starburst galaxy sequence. Second panel: T<sfr 
vs. Y> g /Tdyn- Since IDaddi et al.l (120101 ) report a single relation here, all points from that 
paper have the same style. The Lyman-a galaxies are consistent with the general relation, 
though higher at an insignificant but intriguing level (~ 2a). 
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by lDaddi et all ( 120101 )). This possible deviation should be explored using a larger sample. 



The range of gas surface densities plotted in figure [2] is from ~ 30 to ~ 1 000 M M pc 



-2 



Base d on a standard ratio of dust to gas column density, As ~ S s / (11M pc 2 ) (jBohlin et al. 



19781 ). this corresponds to Ab ~ 3-90 magnitudes of extinction. Yet, our SED fits suggest 



modest extinctions, A v < 1, in all cases. There are a few possible explanations. First, the 
gas surface density could be much lower than one would expect for the observed level of 
star formation activity. In this case, the star formation vs. gas surface density scaling must 
be more extreme than even the starburst relation. Second, the dust-to-gas ratio could be 
about 1-2 orders of magnitude lower than in the Milky Way. This would be most easily 
accommodated if the dust-to-gas ratio scales as the square of metal abundance. Third, 
the extinctions inferred from SED fitting could be dramatic under-estimates. If so, these 
objects' bolometric luminosity would mostly emerge in the rest-frame far infra red, making 
them readily detectible with submillimeter imaging (cf. iFinkelstein et al.l ( 120091 )). 



Densities 



Given our estimates of the dynamical mass, it is straightforward to determine the mean 
density within the effective radius for our sample. We find 



a 



2nG rl 



3vr 

Gt lrb 



where we have used v c = \r2cr and t or b = 2nr e /v c . 

Hierarchical structure formation models suggest that most collapsed galaxies at z ~ 3 
should be incorporated into early-type galaxies or disk galaxy bulges by the present epoch. 
We therefore compare the density measurements for Lyman-a galaxies to corresponding 
densities for early-type galaxies and bulges in figure [3j The Lyman-a galaxy densities are 
the mean density within the effective radius, derived dynamically and therefore inclusive of 
all gravitating mass. Fo r disk galaxy bulges, we take an intermediate-redshift sample from 
MacArthur et al.l (120081 ) . Here we use the published effective radii and velocity dispersions 
to ob tain dynamic a l estim ates of density. We use two nearby elliptical g alaxy samples (one 
from Kelson et al. ( hoOO ). and one compiled by Bezanson et al. ( 2009 ) from earlier work 



bvlFranx et al.l (Il989[ ): IPeletier et al.l ( 1l990l ): l.ledrzejewskil ( 119871 )). For the Kelson et al. 



(2000) sample, we det ermine densities from published a and r e measurements; while for the 
Bezanson et al.l ( 120091 ) sample, we use plotted values of r e and p{r e ). The Lyman-a galaxies 
are smaller at fixed density, and less dense at fixed radius, than are the nearby ellipticals 
and spiral bulges. On the other hand they tend to be smaller and denser than the handful 
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of local ellipticals that have si milar dynamical mas ses. We also compare to a sample of 
early type galaxies at z ~ 2.3 (IBezanson et al.l 120091 ). which overlaps our Lyman-a sample 
in re dshift. Here, the Lyma n-q galaxies are of considerably lower density. In the samples 
from IBezanson et al.l ( 120091 ). we are using stellar mass rather than dynamical mass, since 
a is generally unavailable for the high-redshift early type galaxies. Were dynamical masses 
available, we might expect these two samples to shift upwards in figure [3j but the effect 
would be modest, since the central regions of early type galaxies are likely dominated by 
mass in stars. 

The relation between initial overdensity and epoch of collapse ensures that bound galax- 
ies at 2 < z < 3 will, by the present epoch, be parts of more massive structures. We expect 
that the Lyman-a galaxies we observe will grow through some combination of merging and 
smooth accretion. Merging can proceed either without dissipation, as expected for pure 
mergers of stellar systems ("dry mergers"); or with dissipation, as generally expected for 
gas-rich objects. Conservation of energy and the virial theorem can be combined to infer 
the expected evolution of radius with mass for dry mergers, both in the limit of equal-mass 
("major") mergers, and in the limit of large mass ratio ("minor") mergers. For major m erg- 



ers, we expect R e oc M, while for minor, we expect R e oc M 2 (e.g., IBezanson et al.ll2009l ). In 



figure [31 these correspond to p oc R~ 2 for major mergers, and p oc R~ 2 ' 5 for minor. We plot 
vectors for both scalings. These show that the Lyman-a galaxies we observe cannot grow 
to reproduce the observed masses, sizes, and densities of modern elliptical galaxies through 
simple dry merging. Instead, dissipational ( "wet" ) merging is required, in order to grow the 
galaxies at approximately constant mass density. 



8. Conclusions 



We have performed the first study of dynamical masses for h igh redshift Lyman-o ; 



galaxies, using a sample of nine objects assembled from four samples (iMcLinden et al. 



2012; 



Nilsson et al. 



Richardson et all2013. in prepl ; iFinkelstein et aU201ll ; lHashimoto et al 



2012 



2011 



Nakajima et al 



201lh . Such studies have been previously impractical , given small sizes 



that require HST imaging for size measurements (IMalhotra et al.ll2012l ). faint continuum 
that precludes absorption line measurements, and redshifts that require near-infrared spec- 
troscopy to access those emission lines most useful for kinematic line width measurements. 
We measure dynamical masses ranging from 1O 9 M to 1O 1O M . 

We combine these dynamical masses with stellar masses, to study the position of Lyman- 
a galaxies on a version of the stellar mass Tully-Fisher relation between mass and line width. 
To derive stellar masses consistently for the full sample, we fitted population models to 
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Fig. 3. — Densities and effective radii for Lyman-a galaxies (green circl es) and compari 



son sa mples. Cyan triangles are intermediate-redshift disk gala xy bulges iMacArthur et al. 



( 2008 ). and red squares are low-redshift elliptical galaxies from Kelson et al. ( 200ct ). both 
with densities derived dynamically. Blue asterisks are high-redshift (z ~ 2.3) early-type 
galaxies, and t h ree-p ointed stars are another local early type galaxy sample, both from 
Bezanson et all ( 120091 ). and both using stellar masses estimates from SED fitting. Diagonal 
dotted lines mark masses of 10 9 , 10 10 , 10 11 , and 10 12 M Q , and dashed lines mark surface 
densities of 10 8 , 10 9 , and 10 10 M Q kpc~ 2 . The Lyman-a galaxies are typically about lOx less 
massive than the local ellipticals, yet of comparable density. Non-dissipative ( "dry" ) merging 
results in galaxy densities that decrease as p oc R k e fo r 2 < k < 2.5, based on virial theorem 
arguments and conservation of total energy (see, e.g., iBezanson et all (120091 )). Arrows show 
such scalings, with the shorter arrow corresponding to one major merger and the longer to a 
mass doubling through minor mergers. Neither arrow approaches the region of the elliptical 
galaxy samples, indicating that dissipational merging is required if these objects are indeed 
progenitors of present day early type galaxies. 
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extensive multiband photometry covering rest wavelengths from Lyman-a to beyond the 
4000 A break. This affords constraints on young stars, older stars, and dust. 

We find that the Lyman-a galaxies are broadly consistent with the stellar mass Tully- 
Fisher relation established at lower redshift. Thus, whatever physical conditions allow the 
production and escape of significant Lyman-a, they do not result in strong departures from 
the linewidth-mass relation. On the other hand, about half the sample galaxies are consistent 
with stellar masses significantly below their dynamical masses. In these cases, there is a 
possibility of a dynamically significant reservoir of gas that is present in the inner kpc region 
and provides fuel for ongoing star formation. The dynamical mass measurements are in fact 
the most powerful constraint we have on the presence of old stars in the sample galaxies. 
In most cases, the stellar population fits allow up to ~ 3 x 10 10 M Q of old stars before 
the resulting light appreciably degrades the quality of the SED fit. However, dynamical 
constraints place a tighter limit of < 1 x 1O 1O M on the total mass for the galaxies in our 
sample. 

By using the dynamical masses as an upper bound on the gas mass, and using the 
Ha line measurements and/or stellar population fits to infer star formation rates, we have 
examined the way the Lyman-a galaxies fall on the scaling relations for star formation. 
We conclude that they form stars more actively than the "normal" star forming galaxy 
population at comparable gas mass surface density. Their behavior is consistent with that 
observed in starburst galaxies, despite the typically smaller masses and sizes of the Lyman-a 
galaxy population. 

The dynamical masses we infer remain 1-2 orders of magnitude below the characteristic 
dark halo masses inferred from the clustering of Lyman-a galaxies. Since the ratio of baryonic 
to dark mass should be globally uniform at around 14%, we infer that a significant part of 
the baryonic matter associated with these halos has not yet accreted to the central kpc, 
where the active star formation is observed. 

Finally, we examined densities of these objects. While hierarchical structure formation 
models suggest that small galaxies at z ~ 2-3 should become parts of early type galaxies 
or spiral bulges by z = 0, the Lyman-a galaxies are of smaller size and comparable density 
to present day elliptical galaxies, and of comparable size but smaller density than present 
day spiral galaxy bulges. If the Lyman-a galaxies are to evolve into either, they must do so 
through dissipational merging. This is consistent with the picture of Lyman-a galaxies as 
young, starbursting galaxies, whose present properties and future evolution include a large 
role for gas physics. 
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Table 1: Basic observational properties of the sample. 



ID 


z 


a 


r e (kpc) 


SFR 


Refs.^ 


LAE40844 


3.1117 


81 


1.1 


120 


1 


LAE27878 


3.11879 


43 


1.3 


34 


1 


HPS194 


2.287 


64.7 


1.5 


18 


2 


HPS256 


2.491 


74.0 


1.1 


18 


2 


COSMOS13636 


2.1621 


57 


0.787 


38 


3 


COSMOS30679 


2.19855 


55 


1.825 


395 


3 


CDFS3865 


2.173 


105 


0.955 


178 


3 


CDFS6482 


2.205 


42 


1.78 


71 


3 


LAE-COSMOS-47 


2.24654 


30 


1.139 


9.2 


4 



a References — 1. McLinden et al. (2011) ; Richardson et al. f2013. in prep) ; 2. 
Finkelstein et al. (2011); Song fc al (2013. in prep ): 3. Hashimoto et al. (2012) ; 
Nakaiima et al. (2012); 4. Nilsson et al. (2011) 



Table 2. Best fitting SED model for each object. 



ID Age l a Mass 1 Age 2 a Mass 2 Age 3 a Mass 3 M tot A B \ 2 N,, ands b 



LAE40844 


2 


2.78 x 10 8 


8 


5.53 x 10 s 


300 


5.83 x 10 9 


6.67 x 10 9 


0.217 


42.227 


11 


LAE27878 


2 


6.96 x 10 7 


600 


9.96 x 10 9 








1.003 x 10 10 


0.164 


3.269 


11 


HPS194 


2 


5.52 x 10 7 


300 


2.18 x 10 10 








2.185 x 10 10 





14.92 


9 


HPS256 


2 


3.82 x 10 7 


1500 


3.81 x 10 9 








3.84 x 10 9 


0.110 


2.656 


9 


COSMOS13636 


4 


1.06 x 10 s 


300 


8.53 x 10 9 








8.63 x 10 9 


0.258 


35.15 


9 


COSMOS30679 


4 


1.00 x 10 9 


150 


3.95 x 10 9 








4.95 x 10 9 


0.906 


20.22 


9 


CDFS3865 


2 


3.24 x 10 s 


150 


7.88 x 10 9 








8.21 x 10 9 


0.211 


7.672 


12 


CDFS6482 


2 


1.44 x 10 8 


16 


6.23 x 10 s 


2000 


2.25 x 10 10 


2.325 x 10 10 


0.299 


3.862 


11 


LAE-COSMOS-47 


4 


2.63 x 10 7 


300 


1.03 x 10 9 








1.20 x 10 9 





27.5 


19 



a Ages in Myr. 

b Number of photometric bands used in SED fitting. Wavelength coverage always extends at least from B band (0.44^mi) to IRAC 
Channel 1 (3.6/xm). 
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Table 3. Range of SED fitting results for the sample. 



ID 


M 






M» 




1(7 


M», 




M, 




2cT 


A B ,-2<r 


A B ,-la 


A BMat 




A B , +2a 


LAE40844 


1.00 


X 


10 9 


2.90 


X 


10 9 


3.62 > 


< 10 10 


4.16 


X 


10 1() 


0.081 


0.150 


0.217 


0.22 


0.241 


LAE27878 


1.77 


X 


10 9 


4.05 


X 


10 9 


2.39 x 10 10 


3.85 


X 


10 1() 








0.164 


0.201 


0.254 


HPS194 


2.05 


X 


10 10 


2.12 


X 


10 10 


2.26 x 10 10 


2.45 


X 


10 10 











0.030 


0.100 


HPS256 


1.03 


X 


10 8 


1.87 


X 


10 9 


1.20 > 


< 10 10 


2.16 


X 


10 10 





0.009 


0.110 


0.199 


0.298 


COSMOS13636 


2.61 


X 


10 9 


4.49 


X 


10 9 


1.17 > 


< 10 10 


2.20 


X 


10 10 


0.083 


0.161 


0.258 


0.291 


0.320 


COSMOS30679 


1.00 


X 


10 9 


1.75 


X 


10 9 


1.28 > 


< 10 10 


2.82 


X 


10 1() 


0.718 


0.831 


0.906 


0.991 


1.141 


CDFS3865 


6.72 


X 


10 9 


6.94 


X 


10 9 


1.52 x 10 10 


1.74 


X 


10 1() 


0.184 


0.196 


0.211 


0.230 


0.259 


CDFS6482 


2.10 


X 


10 9 


9.68 


X 


10 9 


3.48 x 10 10 


3.75 


X 


10 10 


0.097 


0.235 


0.299 


0.332 


0.358 


LAE-COSMOS-47 


5.77 


X 


10 7 


5.77 


X 


10 8 


2.13 x 10 9 


3.33 


X 


10 9 











0.093 


0.172 



Table 4: Masses and densities of the sample. 

ID y^-dyn log(Pe) M t>min M,. 



LAE40844 


6. 


.7 


X 


10 9 


-22. 


.39 


2. 


.90 


X 


10 9 


3 


.62 


X 


10 10 


LAE27878 


2. 


.2 


X 


10 9 


-23 


.08 


4. 


.05 


X 


10 9 


2 


.39 


X 


10 10 


HPS194 


5. 


.8 


X 


10 9 


-22 


,85 


2. 


12 


X 


10 10 


2 


.26 


X 


10 10 


HPS256 


5. 


.6 


X 


10 9 


-22. 


.47 


1. 


.87 


X 


10 9 


1 


.20 


X 


10 10 


COSMOS13636 


2. 


.4 


X 


10 9 


-22 


,40 


4. 


.49 


X 


10 9 


1 


.17 


X 


10 10 


COSMOS30679 


5. 


1 


X 


10 9 


-23 


,16 


1. 


.75 


X 


10 9 


1 


.28 


X 


10 10 


CDFS3865 


9. 


.8 


X 


10 9 


-22. 


.04 


6. 


.94 


X 


10 9 


1 


.52 


X 


10 10 


CDFS6482 


2. 


.9 


X 


10 9 


-23 


.38 


9. 


.68 


X 


10 9 


3 


.48 


X 


10 10 


LAE-COSMOS-47 


9. 


.5 


X 


10 s 


-23 


.28 


5. 


.77 


X 


10 s 


2. 


,13 


X 


10 9 
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